Selection, optimization and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse US populations.

Polygenic risk scores (PRSs) have improved in predictive performance, but several challenges remain to be addressed before PRSs can be implemented in the clinic, including reduced predictive performance of PRSs in diverse populations, and the interpretation and communication of genetic results to both providers and patients. To address these challenges, the National Human Genome Research Institute-funded Electronic Medical Records and Genomics (eMERGE) Network has developed a framework and pipeline for return of a PRS-based genome-informed risk assessment to 25,000 diverse adults and children as part of a clinical study. From an initial list of 23 conditions, ten were selected for implementation based on PRS performance, medical actionability and potential clinical utility, including cardiometabolic diseases and cancer. Standardized metrics were considered in the selection process, with additional consideration given to strength of evidence in African and Hispanic populations. We then developed a pipeline for clinical PRS implementation (score transfer to a clinical laboratory, validation and verification of score performance), and used genetic ancestry to calibrate PRS mean and variance, utilizing genetically diverse data from 13,475 participants of the All of Us Research Program cohort to train and test model parameters. Finally, we created a framework for regulatory compliance and developed a PRS clinical report for return to providers and for inclusion in an additional genome-informed risk assessment. The initial experience from eMERGE can inform the approach needed to implement PRS-based testing in diverse clinical settings.

Protein interaction networks in the vasculature prioritize genes and pathways underlying coronary artery disease.

Population-based association studies have identified many genetic risk loci for coronary artery disease (CAD), but it is often unclear how genes within these loci are linked to CAD. Here, we perform interaction proteomics for 11 CAD-risk genes to map their protein-protein interactions (PPIs) in human vascular cells and elucidate their roles in CAD. The resulting PPI networks contain interactions that are outside of known biology in the vasculature and are enriched for genes involved in immunity-related and arterial-wall-specific mechanisms. Several PPI networks derived from smooth muscle cells are significantly enriched for genetic variants associated with CAD and related vascular phenotypes. Furthermore, the networks identify 61 genes that are found in genetic loci associated with risk of CAD, prioritizing them as the causal candidates within these loci. These findings indicate that the PPI networks we have generated are a rich resource for guiding future research into the molecular pathogenesis of CAD.

Selection, optimization, and validation of ten chronic disease polygenic risk scores for clinical implementation in diverse populations.

Polygenic risk scores (PRS) have improved in predictive performance supporting their use in clinical practice. Reduced predictive performance of PRS in diverse populations can exacerbate existing health disparities. The NHGRI-funded eMERGE Network is returning a PRS-based genome-informed risk assessment to 25,000 diverse adults and children. We assessed PRS performance, medical actionability, and potential clinical utility for 23 conditions. Standardized metrics were considered in the selection process with additional consideration given to strength of evidence in African and Hispanic populations. Ten conditions were selected with a range of high-risk thresholds: atrial fibrillation, breast cancer, chronic kidney disease, coronary heart disease, hypercholesterolemia, prostate cancer, asthma, type 1 diabetes, obesity, and type 2 diabetes. We developed a pipeline for clinical PRS implementation, used genetic ancestry to calibrate PRS mean and variance, created a framework for regulatory compliance, and developed a PRS clinical report. eMERGE's experience informs the infrastructure needed to implement PRS-based implementation in diverse clinical settings.

Using brain cell-type-specific protein interactomes to interpret neurodevelopmental genetic signals in schizophrenia.

Genetics have nominated many schizophrenia risk genes and identified convergent signals between schizophrenia and neurodevelopmental disorders. However, functional interpretation of the nominated genes in the relevant brain cell types is often lacking. We executed interaction proteomics for six schizophrenia risk genes that have also been implicated in neurodevelopment in human induced cortical neurons. The resulting protein network is enriched for common variant risk of schizophrenia in Europeans and East Asians, is down-regulated in layer 5/6 cortical neurons of individuals affected by schizophrenia, and can complement fine-mapping and eQTL data to prioritize additional genes in GWAS loci. A sub-network centered on HCN1 is enriched for common variant risk and contains proteins (HCN4 and AKAP11) enriched for rare protein-truncating mutations in individuals with schizophrenia and bipolar disorder. Our findings showcase brain cell-type-specific interactomes as an organizing framework to facilitate interpretation of genetic and transcriptomic data in schizophrenia and its related disorders.

Genome-wide polygenic score to predict chronic kidney disease across ancestries.

Chronic kidney disease (CKD) is a common complex condition associated with high morbidity and mortality. Polygenic prediction could enhance CKD screening and prevention; however, this approach has not been optimized for ancestrally diverse populations. By combining APOL1 risk genotypes with genome-wide association studies (GWAS) of kidney function, we designed, optimized and validated a genome-wide polygenic score (GPS) for CKD. The new GPS was tested in 15 independent cohorts, including 3 cohorts of European ancestry (n = 97,050), 6 cohorts of African ancestry (n = 14,544), 4 cohorts of Asian ancestry (n = 8,625) and 2 admixed Latinx cohorts (n = 3,625). We demonstrated score transferability with reproducible performance across all tested cohorts. The top 2% of the GPS was associated with nearly threefold increased risk of CKD across ancestries. In African ancestry cohorts, the APOL1 risk genotype and polygenic component of the GPS had additive effects on the risk of CKD.

Development and validation of a trans-ancestry polygenic risk score for type 2 diabetes in diverse populations.

BACKGROUND: Type 2 diabetes (T2D) is a worldwide scourge caused by both genetic and environmental risk factors that disproportionately afflicts communities of color. Leveraging existing large-scale genome-wide association studies (GWAS), polygenic risk scores (PRS) have shown promise to complement established clinical risk factors and intervention paradigms, and improve early diagnosis and prevention of T2D. However, to date, T2D PRS have been most widely developed and validated in individuals of European descent. Comprehensive assessment of T2D PRS in non-European populations is critical for equitable deployment of PRS to clinical practice that benefits global populations. METHODS: We integrated T2D GWAS in European, African, and East Asian populations to construct a trans-ancestry T2D PRS using a newly developed Bayesian polygenic modeling method, and assessed the prediction accuracy of the PRS in the multi-ethnic Electronic Medical Records and Genomics (eMERGE) study (11,945 cases; 57,694 controls), four Black cohorts (5137 cases; 9657 controls), and the Taiwan Biobank (4570 cases; 84,996 controls). We additionally evaluated a post hoc ancestry adjustment method that can express the polygenic risk on the same scale across ancestrally diverse individuals and facilitate the clinical implementation of the PRS in prospective cohorts. RESULTS: The trans-ancestry PRS was significantly associated with T2D status across the ancestral groups examined. The top 2% of the PRS distribution can identify individuals with an approximately 2.5-4.5-fold of increase in T2D risk, which corresponds to the increased risk of T2D for first-degree relatives. The post hoc ancestry adjustment method eliminated major distributional differences in the PRS across ancestries without compromising its predictive performance. CONCLUSIONS: By integrating T2D GWAS from multiple populations, we developed and validated a trans-ancestry PRS, and demonstrated its potential as a meaningful index of risk among diverse patients in clinical settings. Our efforts represent the first step towards the implementation of the T2D PRS into routine healthcare.

Multiancestral polygenic risk score for pediatric asthma.

BACKGROUND: Asthma is the most common chronic condition in children and the third leading cause of hospitalization in pediatrics. The genome-wide association study catalog reports 140 studies with genome-wide significance. A polygenic risk score (PRS) with predictive value across ancestries has not been evaluated for this important trait. OBJECTIVES: This study aimed to train and validate a PRS relying on genetic determinants for asthma to provide predictions for disease occurrence in pediatric cohorts of diverse ancestries. METHODS: This study applied a Bayesian regression framework method using the Trans-National Asthma Genetic Consortium genome-wide association study summary statistics to derive a multiancestral PRS score, used one Electronic Medical Records and Genomics (eMERGE) cohort as a training set, used a second independent eMERGE cohort to validate the score, and used the UK Biobank data to replicate the findings. A phenome-wide association study was performed using the PRS to identify shared genetic etiology with other phenotypes. RESULTS: The multiancestral asthma PRS was associated with asthma in the 2 pediatric validation datasets. Overall, the multiancestral asthma PRS has an area under the curve (AUC) of 0.70 (95% CI, 0.69-0.72) in the pediatric validation 1 and AUC of 0.66 (0.65-0.66) in the pediatric validation 2 datasets. We found significant discrimination across pediatric subcohorts of European (AUC, 95% CI, 0.60 and 0.66), African (AUC, 95% CI, 0.61 and 0.66), admixed American (AUC, 0.64 and 0.70), Southeast Asian (AUC, 0.65), and East Asian (AUC, 0.73) ancestry. Pediatric participants with the top 5% PRS had 2.80 to 5.82 increased odds of asthma compared to the bottom 5% across the training, validation 1, and validation 2 cohorts when adjusted for ancestry. Phenome-wide association study analysis confirmed the strong association of the identified PRS with asthma (odds ratio, 2.71, PFDR = 3.71 × 10-65) and related phenotypes. CONCLUSIONS: A multiancestral PRS for asthma based on Bayesian posterior genomic effect sizes identifies increased odds of pediatric asthma.

Endothelial ARHGEF26 is an angiogenic factor promoting VEGF signalling.

AIMS: Genetic studies have implicated the ARHGEF26 locus in the risk of coronary artery disease (CAD). However, the causal pathways by which DNA variants at the ARHGEF26 locus confer risk for CAD are incompletely understood. We sought to elucidate the mechanism responsible for the enhanced risk of CAD associated with the ARHGEF26 locus. METHODS AND RESULTS: In a conditional analysis of the ARHGEF26 locus, we show that the sentinel CAD-risk signal is significantly associated with various non-lipid vascular phenotypes. In human endothelial cell (EC), ARHGEF26 promotes the angiogenic capacity, and interacts with known angiogenic factors and pathways. Quantitative mass spectrometry showed that one CAD-risk coding variant, rs12493885 (p.Val29Leu), resulted in a gain-of-function ARHGEF26 that enhances proangiogenic signalling and displays enhanced interactions with several proteins partially related to the angiogenic pathway. ARHGEF26 is required for endothelial angiogenesis by promoting macropinocytosis of Vascular Endothelial Growth Factor Receptor 2 (VEGFR2) on cell membrane and is crucial to Vascular Endothelial Growth Factor (VEGF)-dependent murine vessel sprouting ex vivo. In vivo, global or tissue-specific deletion of ARHGEF26 in EC, but not in vascular smooth muscle cells, significantly reduced atherosclerosis in mice, with enhanced plaque stability. CONCLUSIONS: Our results demonstrate that ARHGEF26 is involved in angiogenesis signaling, and that DNA variants within ARHGEF26 that are associated with CAD risk could affect angiogenic processes by potentiating VEGF-dependent angiogenesis.

ZBTB33 is mutated in clonal hematopoiesis and myelodysplastic syndromes and impacts RNA splicing.

Clonal hematopoiesis results from somatic mutations in cancer driver genes in hematopoietic stem cells. We sought to identify novel drivers of clonal expansion using an unbiased analysis of sequencing data from 84,683 persons and identified common mutations in the 5-methylcytosine reader, ZBTB33, as well as in YLPM1, SRCAP, and ZNF318. We also identified these mutations at low frequency in myelodysplastic syndrome patients. Zbtb33 edited mouse hematopoietic stem and progenitor cells exhibited a competitive advantage in vivo and increased genome-wide intron retention. ZBTB33 mutations potentially link DNA methylation and RNA splicing, the two most commonly mutated pathways in clonal hematopoiesis and MDS.

Genoppi is an open-source software for robust and standardized integration of proteomic and genetic data.

Combining genetic and cell-type-specific proteomic datasets can generate biological insights and therapeutic hypotheses, but a technical and statistical framework for such analyses is lacking. Here, we present an open-source computational tool called Genoppi (lagelab.org/genoppi) that enables robust, standardized, and intuitive integration of quantitative proteomic results with genetic data. We use Genoppi to analyze 16 cell-type-specific protein interaction datasets of four proteins (BCL2, TDP-43, MDM2, PTEN) involved in cancer and neurological disease. Through systematic quality control of the data and integration with published protein interactions, we show a general pattern of both cell-type-independent and cell-type-specific interactions across three cancer cell types and one human iPSC-derived neuronal cell type. Furthermore, through the integration of proteomic and genetic datasets in Genoppi, our results suggest that the neuron-specific interactions of these proteins are mediating their genetic involvement in neurodegenerative diseases. Importantly, our analyses suggest that human iPSC-derived neurons are a relevant model system for studying the involvement of BCL2 and TDP-43 in amyotrophic lateral sclerosis.

Targeting acute myeloid leukemia dependency on VCP-mediated DNA repair through a selective second-generation small-molecule inhibitor.

The development and survival of cancer cells require adaptive mechanisms to stress. Such adaptations can confer intrinsic vulnerabilities, enabling the selective targeting of cancer cells. Through a pooled in vivo short hairpin RNA (shRNA) screen, we identified the adenosine triphosphatase associated with diverse cellular activities (AAA-ATPase) valosin-containing protein (VCP) as a top stress-related vulnerability in acute myeloid leukemia (AML). We established that AML was the most responsive disease to chemical inhibition of VCP across a panel of 16 cancer types. The sensitivity to VCP inhibition of human AML cell lines, primary patient samples, and syngeneic and xenograft mouse models of AML was validated using VCP-directed shRNAs, overexpression of a dominant-negative VCP mutant, and chemical inhibition. By combining mass spectrometry-based analysis of the VCP interactome and phospho-signaling studies, we determined that VCP is important for ataxia telangiectasia mutated (ATM) kinase activation and subsequent DNA repair through homologous recombination in AML. A second-generation VCP inhibitor, CB-5339, was then developed and characterized. Efficacy and safety of CB-5339 were validated in multiple AML models, including syngeneic and patient-derived xenograft murine models. We further demonstrated that combining DNA-damaging agents, such as anthracyclines, with CB-5339 treatment synergizes to impair leukemic growth in an MLL-AF9-driven AML murine model. These studies support the clinical testing of CB-5339 as a single agent or in combination with standard-of-care DNA-damaging chemotherapy for the treatment of AML.

Cohesin mutations alter DNA damage repair and chromatin structure and create therapeutic vulnerabilities in MDS/AML.

The cohesin complex plays an essential role in chromosome maintenance and transcriptional regulation. Recurrent somatic mutations in the cohesin complex are frequent genetic drivers in cancer, including myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML). Here, using genetic dependency screens of stromal antigen 2-mutant (STAG2-mutant) AML, we identified DNA damage repair and replication as genetic dependencies in cohesin-mutant cells. We demonstrated increased levels of DNA damage and sensitivity of cohesin-mutant cells to poly(ADP-ribose) polymerase (PARP) inhibition. We developed a mouse model of MDS in which Stag2 mutations arose as clonal secondary lesions in the background of clonal hematopoiesis driven by tet methylcytosine dioxygenase 2 (Tet2) mutations and demonstrated selective depletion of cohesin-mutant cells with PARP inhibition in vivo. Finally, we demonstrated a shift from STAG2- to STAG1-containing cohesin complexes in cohesin-mutant cells, which was associated with longer DNA loop extrusion, more intermixing of chromatin compartments, and increased interaction with PARP and replication protein A complex. Our findings inform the biology and therapeutic opportunities for cohesin-mutant malignancies.

BraInMap Elucidates the Macromolecular Connectivity Landscape of Mammalian Brain.

Connectivity webs mediate the unique biology of the mammalian brain. Yet, while cell circuit maps are increasingly available, knowledge of their underlying molecular networks remains limited. Here, we applied multi-dimensional biochemical fractionation with mass spectrometry and machine learning to survey endogenous macromolecules across the adult mouse brain. We defined a global "interactome" comprising over one thousand multi-protein complexes. These include hundreds of brain-selective assemblies that have distinct physical and functional attributes, show regional and cell-type specificity, and have links to core neurological processes and disorders. Using reciprocal pull-downs and a transgenic model, we validated a putative 28-member RNA-binding protein complex associated with amyotrophic lateral sclerosis, suggesting a coordinated function in alternative splicing in disease progression. This brain interaction map (BraInMap) resource facilitates mechanistic exploration of the unique molecular machinery driving core cellular processes of the central nervous system. It is publicly available and can be explored here https://www.bu.edu/dbin/cnsb/mousebrain/.

Genome-scale analysis identifies paralog lethality as a vulnerability of chromosome 1p loss in cancer.

Functional redundancy shared by paralog genes may afford protection against genetic perturbations, but it can also result in genetic vulnerabilities due to mutual interdependency1-5. Here, we surveyed genome-scale short hairpin RNA and CRISPR screening data on hundreds of cancer cell lines and identified MAGOH and MAGOHB, core members of the splicing-dependent exon junction complex, as top-ranked paralog dependencies6-8. MAGOHB is the top gene dependency in cells with hemizygous MAGOH deletion, a pervasive genetic event that frequently occurs due to chromosome 1p loss. Inhibition of MAGOHB in a MAGOH-deleted context compromises viability by globally perturbing alternative splicing and RNA surveillance. Dependency on IPO13, an importin-ß receptor that mediates nuclear import of the MAGOH/B-Y14 heterodimer9, is highly correlated with dependency on both MAGOH and MAGOHB. Both MAGOHB and IPO13 represent dependencies in murine xenografts with hemizygous MAGOH deletion. Our results identify MAGOH and MAGOHB as reciprocal paralog dependencies across cancer types and suggest a rationale for targeting the MAGOHB-IPO13 axis in cancers with chromosome 1p deletion.

Water Freezing and Ice Melting.

The generalized replica exchange method (gREM) is designed to sample states with coexisting phases and thereby to describe strong first order phase transitions. The isobaric MD version of the gREM is presented and applied to the freezing of liquid water and the melting of hexagonal and cubic ice. It is confirmed that coexisting states are well-sampled. The statistical temperature as a function of enthalpy, TS(H), is obtained. Hysteresis between freezing and melting is observed and discussed. The entropic analysis of phase transitions is applied and equilibrium transition temperatures, latent heats, and surface tensions are obtained for hexagonal ice ↔ liquid and cubic ice ↔ liquid with excellent agreement with published values. A new method is given to assign water molecules among various symmetry types. Pathways for water freezing, ultimately leading to hexagonal ice, are found to contain intermediate layered structures built from hexagonal and cubic ice.

Pathways through Equilibrated States with Coexisting Phases for Gas Hydrate Formation.

Under ambient conditions, water freezes to either hexagonal ice or a hexagonal/cubic composite ice. The presence of hydrophobic guest molecules introduces a competing pathway: gas hydrate formation, with the guests in clathrate cages. Here, the pathways of the phase transitions are sought as sequences of states with coexisting phases, using a generalized replica exchange algorithm designed to sample them in equilibrium, avoiding nonequilibrium processes. For a dilute solution of methane in water under 200 atm, initializing the simulation with the full set of replicas leads to methane trapped in hexagonal/cubic ice, while gradually adding replicas with decreasing enthalpy produces the initial steps of hydrate growth. Once a small amount of hydrate is formed, water rearranges to form empty cages, eventually transforming the remainder of the system to metastable ß ice, a scaffolding for hydrates. It is suggested that configurations with empty cages are reaction intermediates in hydrate formation when more guest molecules are available. Free energy profiles show that methane acts as a catalyst reducing the barrier for ß ice versus hexagonal/cubic ice formation.

Isobaric Molecular Dynamics Version of the Generalized Replica Exchange Method (gREM): Liquid-Vapor Equilibrium.

A prescription for sampling isobaric generalized ensembles with molecular dynamics is presented and applied to the generalized replica exchange method (gREM), which was designed to simulate first-order phase transitions. The properties of the isobaric gREM ensemble are discussed, and a study is presented for the liquid-vapor equilibrium of the guest molecules given for gas hydrate formation with the mW water model. Phase diagrams, critical parameters, and a law of corresponding states are obtained.

Entropic Description of Gas Hydrate Ice-Liquid Equilibrium via Enhanced Sampling of Coexisting Phases.

Metastable ß ice holds small guest molecules in stable gas hydrates, so its solid-liquid equilibrium is of interest. However, aqueous crystal-liquid transitions are very difficult to simulate. A new molecular dynamics algorithm generates trajectories in a generalized NPT ensemble and equilibrates states of coexisting phases with a selectable enthalpy. With replicas spanning the range between ß ice and liquid water, we find the statistical temperature from the enthalpy histograms and characterize the transition by the entropy, introducing a general computational procedure for first-order transitions.

Role of Charge and Solvation in the Structure and Dynamics of Alanine-Rich Peptide AKA2 in AOT Reverse Micelles.

The propensity of peptides to form α-helices has been intensely studied using theory, computation, and experiment. Important model peptides for the study of the coil-to-helix transition have been alanine-lysine (AKA) peptides in which the lysine residues are placed on opposite sides of the helix avoiding charge repulsion while enhancing solubility. In this study, the effects of capped versus zwitterionic peptide termini on the secondary structure of alanine-rich peptides in reverse micelles are explored. The reverse micelles are found to undergo substantial shape fluctuations, a property observed in previous studies of AOT reverse micelles in the absence of solvated peptide. The peptides are observed to interact with water, as well as the AOT surfactant, including interactions between the nonpolar residues and the aliphatic surfactant tails. Computation of IR spectra for the amide I band of the peptide allows for direct comparison with experimental spectra. The results demonstrate that capped AKA2 peptides form more stable α helices than zwitterionic AKA2 peptides in reverse micelles. The rotational anisotropy decay of water is found to be distinctly different in the presence or absence of peptide within the reverse micelle, suggesting that the introduction of peptide significantly alters the number of free waters within the reverse micelle nanopool. However, neither the nature of the peptide termini (capped or charged) nor the degree of peptide helicity is found to significantly alter the balance of interactions between the peptides and the environment. Observed changes in the degree of helicity in AKA2 peptides in bulk solution and in reverse micelle environments result from changes in peptide confinement and hydration as well as direct nonpolar and polar interactions with the water-surfactant interface.